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You need not be afraid of phantom energy 
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Phantom energy which violates the dominant-energy condition and is not excluded by current 
constraints on the equation of state may be dominating the evolution of the universe now. It has 
been pointed out that in such a case the fate of the universe may be a big rid where the expansion is 
so violent that all galaxies, planet and even atomic nuclei will be successively ripped apart in finite 
time. Here we show however that there are certain unified models for dark energy which are stable 
to perturbations in matter density where the presence of phantom energy does not lead to such a 
cosmic doomsday. 

PACS numbers: 98.80.-k, 98.80.Es 



I. INTRODUCTION 

WMAP [1] has confirmed it with the highest accuracy: 
Nearly seventy percent of the energy in the universe is in 
the form of dark energy -possibly one of most astonishing 
discoveries ever made in science. Moreover, recent obser- 
vations do not exclude, but actually suggest a value even 
smaller than -1 for the parameter of the equation of state, 
ui, characterizing that dark energy [2]. That means that 
for at least a perfect-fluid equation of state, the absolute 
value for negative pressure exceeds that for positive en- 
ergy density, i.e. p + p < 0, and hence it follows that 
the involved violation of the dominant-energy condition 
might allow the existence of astrophysical or cosmologi- 
cal wormholes. A most striking consequence from dark 
energy with u < — 1 has been however pointed out [3]. It 
is that in a finite time the universe will undergo a catas- 
trophic "big rip" . Big rip is a term coined by Caldwell [3] 
that corresponds to a new cosmological model in which 
the scale factor blows up in a finite time because its cos- 
mic acceleration is larger than that is induced by a cos- 
mological constant, making in this way every component 
of the universe to go beyond the horizon of all other uni- 
verse's components in finite cosmic time. If dark energy 
is phantom energy, i.e. if dark energy is characterized by 
an equation of state with u = p/ p < — 1, and hence there 
is a violation of the dominant-energy condition, p+p < 0, 
then the phantom-energy density is still positive though 
it will first increase from a finite small initial value up 
to infinite in a finite time, thereafter steadily decreas- 
ing down to zero as time goes to infinity. A state with 
infinite energy density at finite cosmological time is cer- 
tainly an unusual state in cosmology. To an observer on 
the Earth, this state coincides with the above-mentioned 
big rip where the universe dies after ripping successively 
apart all galaxies, our solar system, the Earth itself, and 
finally molecules, atoms, nuclei and nucleons [3]. For 
a general cosmological model with phantom energy, the 
time at which that big rip would take place depends on 
both the initial size of the universe and the value of u in 
such a way that, the larger the absolute value of these 
quantities, the nearer the big rip will occur. The be- 



haviour of the universe after the big rip is in some re- 
spects even more bizarre than the big rip itself, as its 
size then steadily decreases from infinite down to zero at 
infinite time. In case that a generic perfect-fluid equation 
of state, p — up, with u < — 1 is considered, the above 
new behaviours show themselves immediately. For flat 
geometry, the scale factor is then given by [4]: 
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with C\ > and < C 2 < 1. We note that for u < -1, 
in fact a — > oo as t — > — lnC?2/[Ci(l + This marks 
the time at big rip and the onset of the contracting phase 
for t > t*. 

This report aims at showing that the above-described 
emergence of a cosmic doomsday at which the big rip 
occurred and the subsequent unconventional evolution 
of the universe can both be avoided while keeping the 
phantom energy condition, p + p <0, p > 0, u < — 1, on 
the dark energy if, instead of a quintessential description 
of dark energy based on an equation of state p = up, 
with u < —1, we consider a suitable generalization of 
the Chaplygin-gas model which, at sufficiently late times, 
does not show observable non-physical oscillations and 
exponential blowup in the matter density perturbations 
[5] that are present in current unstable Chaplygin-gas 
cosmic models [6]. The latter models describe a single 
substance which is characterized by an equation of state 
[6] p = —Ap~ a , where A is a positive-definite constant 
and a is a parameter which may take on any real positive 
values. This equation of state has been shown to rep- 
resent the stiff that simultaneously describes dark mat- 
ter and dark energy, but gives rise to instabilities stem- 
ming from the unobserved existence of oscillations and 
exponential blowup in the perturbation power-spectrum 
which arise whenever the speed of sound in nonzero [5]. 

The paper can be outlined as follows. In Sec. II we 
generalize the cosmic Chaplygin-gas models in such a way 
that the resulting models can be made stable and free 
from unphysical behaviours even when the vacuum fluid 
satisfies the phantom energy condition. The Friedmann 
equations for models which show and do not show un- 
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physical behaviours are solved in Sec. Ill, checking that 
in the latter case the phantom energy condition does not 
imply any emergence of a big rip in finite time. We finally 
conclude in Sec. IV. 



II. GENERALIZED COSMIC CHAPLYGIN-GAS 
MODELS 

We introduce here some generalizations from cosmic 
Chaplygin-gas model that also contain an adjustable ini- 
tial parameter u>. In particular, we shall consider a gen- 
eralized gas whose equation of state reduces to that of 
current Chaplygin unified models for dark matter and 
energy in the limit lu — > and satisfies the following con- 
ditions: (i) it becomes a de Sitter fluid at late time and 
when u) = — 1, (ii) it reduces to p — top in the limit that 
the Chaplygin parameter A — > 0, (iii) it also reduces to 
the equation of state of current Chaplygin unified dark 
matter models at high energy density, and (iv) the evo- 
lution of density perturbations derived from the chosen 
equation of state becomes free from the above-mentioned 
pathological behaviour of the matter power spectrum for 
physically reasonable values of the involved parameters, 
at late time. We shall see that these generalizations re- 
tain a big rip if they also show unphysical oscillations and 
exponential blowup leading to instability (i.e. if they do 
not satisfy condition (iv)), but if the latter effects are 
avoided then the evolution of the scale factor recovers a 
rather conventional pattern, without any big rip or con- 
tracting phase. 

An equation of state that can be shown to satisfy all 
the above conditions (i) - (iv) is 
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with A a constant which now can take on both positive 
and negative values, and > u > —£, I being a positive 
definite constant which can take on values larger than 
unity. By integrating the cosmic conservation law for 
energy we get for the energy density 
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where B is a positive integration constant. Let us now 
define the effective expressions of the state equation pa- 
rameter and speed of sound, which respectively are given 
by 
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One can then interpret the model by taking the limit 
of these parameters as a — > and a — > oo, at which 
limits they respectively become iv cfi — > and c° s S2 — > 
(that correspond to the pressureless CDM model), and 
iv cfi — + — 1 and c c s S2 — > a + u(l + a) (that correspond to a 
pure cosmological constant). Such as we have defined it 
so far, the present model does not satisfy condition (iv) 
above, as the evolution of density perturbations 6 k with 
wave vector k [5] 

6' k ' + F(u; cS , cf)8' k G(c off , cf, k)5 k = 0, (2.7) 

where the prime denotes differentiation with respect to 
In a, and 
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with H the Hubble parameter, 
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and 1]m the CDM density defined from Eq. (2.3) in the 
limit a — > 0, shows oscillations and exponential blowup 
because c° s is generally nonzero. 



III. AVOIDING THE BIG RIP 

In case that u> < — 1 for the large values of the scale 
factor for which dominance of the dark-phantom energy 
is expected, one can approximate the Friedmann equa- 
tion that corresponds to the considered model as follows 
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in which I? = STrGB-Via+aXM-i)]. The solution to 
this equation is 
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where ao and to are the initial values of the scale factor 
and time, respectively. It is easy to see that in the consid- 
ered case the phantom energy condition always satisfies 
p + p < and there will be a big rip, taking place now 
at a time 
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followed as well by a contracting phase where the size of 
the universe vanishes as t — > oo. Note that, as it was 
pointed out before, the time at which the big rip occurs 
turns out to depend on the initial size of the universe ao, 
in such a way that the big rip becomes nearer as ao is 
made larger. Thus, the above model, which is actually 
excluded because it shares the same kind of instabilities 
as the original cosmic Chaplygin-gas model [5], shows a 
cosmological big rip. 

In order to allow for both stability and compatibility 
with observations, we consider next a model in which 
c° ff — > as t — > oo and the nonzero value of parameter 
B is small enough. The first of these conditions can be 
achieved by simply imposing a + + a) = 0, i.e. 
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The equation of state and the expression for the energy 
density are then reduced to read: 
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We notice that a — > ao as t — > to and a — > oo as t 
and hence there is not a big rip for this solution. 
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This equation of state satisfies then all the condition (i) 
- (iv) imposed above. 

According to condition (3.4) the phantom-dark energy 
with u> < — 1 immediately implies that a < — 1 too. In 
such a case we can check that p + p < for all values 
of the scale factor and, in order to ensure positiveness 
of the energy density, we must have A = —\A\ and keep 
B > 0. Then, for the large values of the scale factor for 
which phantom-dark energy is expected to dominate over 
matter, the energy density can be approximated to 
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We can now write for the Friedmann equation 
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IV. CONCLUSION 

It appears then that if we choose a general equation 
of state for dark energy which is reasonably free from in- 
stabilities and unphysical effects, then a phantom energy 
can be predicted which does not show any big rip at finite 
time. The key difference between the scale factor given 
by Eq. (3.11) and that given by Eq. (1.1) is in the sign 
of the overall exponent of the right-hand-side; while in 
Eq. (1.1) it is negative for uo < —1, in Eq. (3.11) it is 
positive for the same case. Thus, cosmology can co-exist 
with these phantoms in a quite safe manner. 
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